Monte Carlo studies of Matrix theory correlation functions 
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We study correlation functions in (0 + l)-dimensional maximally supersymmetric U(7V) gauge 
theory, which represents the low-energy effective theory of DO-branes. In the large-TV limit, the 
gauge-gravity duality predicts power-law behaviors in the infrared region for the two-point corre- 
lation functions of operators corresponding to supergravity modes. We evaluate such correlation 
functions on the gauge theory side by the Monte Carlo method. Clear power-law behaviors are ob- 
served at TV = 3, and the predicted exponents are confirmed consistently. Our results suggest that 
the agreement extends to the M-theory regime, where the supergravity analysis in 10 dimensions 
may not be justified a priori. 

PACS numbers: ll.25.-w; 11.25.Sq 

Introduction. — Maximally supersymmetric Yang- 
Mills theories (SYM) in various dimensions are impor- 
tant because of their connection to non-perturbative for- 
mulations of string theory. Of particular interest is the 
(0+l)D SYM with U(N) gauge symmetry, which is sup- 
posed to describe the low-energy gravitational dynamics 
of DO-branes in 10D type IIA superstring theory. In an 
appropriate large- N limit, this theory was proposed to 
be a definition of M-theory in a special light-like frame, 
and it is commonly referred to as the Matrix theory [l|. 
M-theory is a hypothetical 11D theory [2j, whose low en- 
ergy effective theory is given by 1 ID supergravity, and it 
is believed to appear in the strong coupling limit of 10D 
type IIA superstring theory. 

Indeed it was confirmed that scattering amplitudes in 
the Matrix theory at weak coupling are consistent with 
predictions from 1 ID supergravity. Even the three-body 
force, which is a characteristic non-linear effect of general 
relativity, has been reproduced 0] from the quantum loop 
effects of massive open strings connecting DO-branes. On 
the other hand, direct perturbative calculations of corre- 
lation functions, which could provide crucial information 
on the as yet mysterious theory, are plagued by severe in- 
frared divergences caused by the massless modes inherent 
in the theory. We therefore need genuinely nonperturba- 
tive methods to study such quantities. 

In the case of A/*=4 SYM in (3+1) dimensions, various 
useful insights have been gained from the AdS/CFT cor- 
respondence 0, HI- This conjectural duality enables us 
to study the M = 4 SYM in the large- N limit with large 
't Hooft coupling constant by using the weakly coupled 
supergravity, which describes the low-energy limit of the 
string theory. Likewise, SYM in (p+l)-dimensions, cor- 
responding to the world-volume theory of Dp-branes, is 
expected to be dual to a superstring theory on the Dp- 



brane background in the near-horizon limit 0, @] • 

Assuming this duality for the p = case and using 
the known dictionary [8(] between the supergravity modes 
and the Matrix theory operators, one may hope to calcu- 
late the correlation functions in the large- N limit, on the 
basis of the Gubser-Klebanov-Polyakov-Witten prescrip- 
tion This program was carried out a decade ago by 
Y.S. and T.Y. Q using the "generalized conformal sym- 
metry" 0, US] as a guide for classifying and interpreting 
the obtained results. The prediction is that the two-point 
correlation functions obey the power law 
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for a set of operators 0(t) corresponding to the super- 
gravity modes. The exponents v are fractional numbers 
which differ from the canonical values, and they are re- 
lated to the generalized conformal dimensions A of the 
operators under consideration byA = + Q. The 
power-law behavior conforms to the existence of a unique 
threshold bound state at zero energy [ll| in the gauge 
theory, which corresponds to single-body states in the 
graviton supermultiplet. 
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This result has been re-derived in later works 
tending the analysis in Ref. Q to general Dp-branes (p < 
5) in the PP-wave limit. There it has also been shown 
that the infrared behavior of the correlation functions for 
operators corresponding to excited stringy modes is typi- 
cally exp{— c(gy M N\t — i'l 3 ) 1 ^ 5 } up to power corrections. 

The near-horizon limit of the DO-brane background is 
valid when the radial coordinate r in the bulk supergrav- 
ity satisfies r -C (gsN) 1 ^ 7 in the string unit a' = 1. The 
analysis using the supergravity approximation is then 
possible when g s N ^> 1, and is trustable when both 
the background curvature and the effective string cou- 
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pling given by the dilaton expectation value are small. 
The string coupling becomes stronger towards the center, 
while the curvature becomes larger towards the bound- 
ary. Taking into account the fact that the radial coordi- 
nate r is related to the time scale t at the boundary as 
r ~ (i/A 1/2 )" 2/5 @] with A = g\ M N ~ g s N 7 the con- 
ditions for justifying the supergravity result (P) can be 
summarized as [a, 0, Q 



A" 1 / 3 « \t - t'\ « A~ 1/3 7V 10/2 



(2) 



Hence, the application of the gauge-gravity correspon- 
dence is legitimate when one studies the region which is 
sufficiently infrared compared to the length scale set by 
A. Since the upper bound in (|2|) diverges at large N, 
the infrared behavior of the correlation functions can be 
predicted reliably from the supergravity analysis in the 
large- N limit. 

In this Letter we evaluate the correlation functions 
on the gauge theory side by the Monte Carlo method 
and compare the results with the predictions from 
the gauge-gravity correspondence. Related Monte Carlo 
analyses have been recently applied by the authors in- 
cluding M.H. and J.N. to reproduce the black hole ther- 
modynamics 3, IH and the Schwarzschild radius [T3| of 
the dual geometry. To the best of our knowledge, this 
is the first attempt to compute the Matrix theory cor- 
relation functions from first principles 18|. Our Monte 



Carlo data are consistent with the supergravity predic- 
tions with reasonable accuracy despite the fact that the 
matrix size is as small as N — 3. Possible interpretations 
shall be given later. 

Simulating Matrix theory. — The action of the Matrix 
theory can be written (in the Euclidean convention) as 

NT 13 r 1 1 

S = Tj Q dtTr{-(D t X l f--[X ll X, j \ 2 

+^aAi - ^tpa(li)oip[Xi,tp^]j , (3) 

where X, {i = 1, • • ■ ,9) and tp a (a = 1, • • • , 16) are NxN 
bosonic and fermionic Hermitian matrices, on which the 
covariant derivative D t acts as D t — dt — i[A, • ] with 
the gauge field A. The 16 x 16 matrices ji satisfy the 
Clifford algebra — 25ij. It is convenient to adopt 

units in which A = 1 without loss of generality. This does 
not contradict the conditions for the bulk theory, since 
the gauge theory has no independent length scale other 

2/3 

than the coupling constant g Y M an d hence the strong 
coupling limit amounts to the IR limit. 

In order to put the system on a computer, we have to 
introduce the UV and IR cutoffs appropriately. The ex- 
tent (3 in the Euclidean time direction represents the IR 
cutoff, which should be sufficiently large in order to see 
the correct infrared properties. Since we are not inter- 
ested in the finite temperature behaviors unlike in Refs. 
14l.ll7j|. we impose periodic boundary conditions for both 



bosonic and fermionic matrices, respecting supersymme- 
try. In order to introduce the UV cutoff, we first fix the 
gauge as A(t) = idiag(ai, • • • , ajy), where — ir < a, < ir, 
and then introduce a Fourier mode cutoff A as Xi(t) — 



and i/} a (t) = J2r, 



4>o 



where 



lu = 2it / (3. Integration over the fermionic matrices yields 
a Pfaffian, which could be complex in general. In this 
work, we simply neglect its phase following the argument 
in Rcf. [17] . While we do not have purely theoretical jus- 
tification of our assumption at this moment, the previous 
works [LU IT?] ] provide empirical evidence that the phase 
does not affect the results in the continuum limit. In- 
cluding its effect is possible in principle by the so-called 
reweighting method, but it is not feasible for our purpose, 
since the demand on the statistics would increase expo- 
nentially with the system size due to huge cancellation. 
This is generally called the "sign problem" . The system 
of finite degrees of freedom obtained under the above as- 
sumption can be simulated by using the Rational Hybrid 
Monte Carlo algorithm [2(| as described in Ref. 14 1. 

Correlation functions. — First let us study a series of 
operators ... it {I > 1), which is defined by 



(4) 



where Fij = —i[Xi,Xj] and Str represents the sym- 
metrized trace treating as a single unit. The value of 
v in |T]) predicted in this case is v — 21/5 

Since we are dealing with the Fourier modes in our 
simulation, the two-point functions which are directly ac- 
cessible are those in the momentum space (o[p)0{— p)\, 

where p — ncu. Correlation functions in the real space 
can be obtained by the inverse Fourier transformation 

(O(i)O(O)) = | J2 P (0(P)0(-P)) e lp *- However, the re- 
sults obtained in this way oscillates with the frequency 
of O(Aw) as a result of the Gibbs phenomenon. For 
the series of operators (@|, we find that the momentum- 
space correlator falls off as p~ 2 at large p. Assuming that 
this behavior continues to infinite p, we extend it up to 
p = 1000 lu, and then make the inverse Fourier trans- 
formation, which removes the Gibbs phenomenon com- 
pletely [131 ] . In Fig.[T] we show the real-space correlation 
functions obtained in this way. The predicted power-law 
behavior can be seen in the region 0.5 < t < 1.5, which is 
somewhat wider than the criterion ^ for justifying the 
supergravity analysis. Note that the IR bound of @ is 
\t — t'\ <C 1-69 for N = 3. Considering that 1.5 is close 
to (3/2 — 2, we expect that the deviation observed at 
t > 1.5 is rather due to the IR cutoff. 



Let us next study the operator T 2 +l fe defined by 



rp + 'i 



-Str ((DtXJXjXk) 



(5) 



for which v is predicted to be v = 4/5 [9j. Unlike the 
previous case, the two-point function in the Fourier space 
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FIG. 1: The real-space two-point correlation functions 
{J+(t)J+(0)} (I = 1, 2, 3, 4) are plotted for N = 3 and (3 = 4. 
The UV cutoff is A = 16. The straight lines are fits to the 
predicted power-law behavior. 



does not fall off as p~ 2 at large p. This is likely due 
to the presence of a derivative Dt in (|5|). Furthermore, 
we observe considerable dependence on the UV cutoff 
A. We therefore evaluate the two-point function in the 
momentum space at various A and make an extrapolation 
to A = oo at each p. The leading finite A effects are 
consistent with 1/A as one might expect theoretically, 
and we make an extrapolation with A = 6, 8, 12 for < 
p/(2tt) < 0.6 and A = 8, 12, 16 for p/(2ir) > 0.8. Figured 
shows the results obtained in this way for N = 3, (3 = 5. 

On the other hand, the prediction (fTJ) in the Fourier 
space reads |9j 



6{p)0{-p))~f(jp) + \p\ 2v g{p) 



(6) 



at small p, where f(p) and g(p) are analytic functions 
(g(0) 7^ 0), which cannot be determined by the super- 
gravity analysis. The odd powers of p in f(p) and g(p) 
are forbidden by the time reflection invariance. For the 
present operator, considering that v = 4/5, a few leading 
terms at small p are predicted to be of the form 



n{ P )ni-p) 



bp 2 + c\p\ 2 



(7) 



where a, b and c are constants. Note that the polyno- 
mial terms do not affect the long-distance behavior of 
the correlation function. Treating v as a free parameter 
and fitting our data in < p/(2ir) < 1 to 0, we obtain 
v = 0.80 ± 0.03 as in Fig. [1 

As the last example, let us consider the operator 



T, 



2. 'J 



1 

N 



tr(X^) (ijLj) 



(8) 



for which the supergravity analysis predicts v = —3/5 
9]. The two-point function is predicted to behave as 
\t — t'l 1 / 5 in the real space, which is divergent as \t — t'\ — » 
oo. This IR divergence is reminiscent of the free-field 
behavior \t — t'\ 2 albeit with much weaker power. Here 



FIG. 2: The momentum-space two-point correlation function 

(f£(p)f£(-p)\ is plotted for N = 3 and (3 = 5. The dotted 

line represents a fit to the behavior treating v as a free 
parameter in the range < p/(2n) < 1. The best fit is 
obtained for v = 0.80 ± 0.03. 



we examine the behavior ([6]) for small p in the Fourier 
space. Considering that v = —3/5, a few leading terms 
are predicted to be of the form 



T+ + {p)T+ + {-p) 



a\p\ 2v + b + c\p\ 



2u+2 



(9) 



Note that the leading term is divergent as p — > 0, as 
opposed to the previous cases with vanishing behaviors. 
In Fig. [31 we show a log- log plot of the two-point function. 
We extrapolated our data to A = oo as we did for T+. 
Treating v as a free parameter and fitting our results 
to © in the region 0.3 < p/(2n) < 2, we obtain v = 
—0.61 ± 0.02. The deviation from the behavior ([9]) at 
small p is most likely due to finite IR cutoff effects as one 
can see from the trends with increasing (3. 
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FIG. 3: The momentum-space two-point correlation function 
(t2 + (p)T2 + (-p)) is plotted in the log-log scale for N = 3. 
The dotted line represents a fit to the behavior © treating v 
as a free parameter in the range 0.3 < p/(2n) < 2. The best 
fit is obtained for v = -0.61 ± 0.02. 



Summary and discussion. — In this Letter we pre- 
sented the first Monte Carlo evaluation of two-point cor- 
relation functions of the supergravity operators in the 
Matrix theory. In particular, we observed that the power 
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law predicted by the gauge-gravity correspondence con- 
tinues to hold in the infrared region beyond the naive 
criterion for the validity of the supergravity analysis. 

From the gauge theory point of view, the existence of a 
threshold bound state at zero energy requires the power 
law at sufficiently large \t — 1'\ for any TV (> 2). In fact, 
the existence of a unique threshold bound state for each 
N (> 2) has been argued, based on the calculation of 
the Witten index [2l|. Let us also consider the spectral 
representation of the two-point function 



0(p)6(- P ) 



dnp{n) 



1 



p z + fi z 



(10) 



The power law of the two-point function at small p im- 
plies that the spectral density p{p) behaves as p(/i) ~ 
p? v+1 at small p. The fact that the exponent v increases 
in general [9j with the angular momentum carried by the 
operator is consistent with the continuous mass spectrum 
of the intermediate many-body states composed of the 
zero-energy bound states. 

On the gravity side, let us emphasize that the upper 
bound in (|2|) is merely a sufficient condition for justifying 
the supergravity analysis. Here we recall that the wave 
function of the supergravity fluctuations obtained in Ref. 
Q is exponentially small in the central region. Therefore, 
it is possible that the 10D supergravity approximation 
in this region is not affected by the increasingly large 
effective string coupling at small r and hence at large 
\t — t'\. On the other hand, in the short-distance regime 
\t— 1'\ <C A -1 / 3 , where the curvature becomes large on the 
gravity side, we expect large deviations from the power- 
law behavior |T]). This is clear, in particular, for the cases 
with 2v + 1 > since there should be no UV singularity 
in one dimension. Such deviations are indeed observed 



in the stringy PP-wave analysis [12] on the bulk side. 

Another significant aspect of our results is that the ex- 
ponents obtained from the Monte Carlo data at N = 3 
essentially coincide with the predictions in the large-A^ 
limit. We have also studied the N = 2 case, but the 
exponents are unaltered within numerical uncertainties. 
This suggests that the exponents are independent of N 
without 0(l/iV 2 ) corrections, possibly due to the BPS 
nature of the supergravity modes. Furthermore, since 
the exponents cannot depend on the coupling constant 
A, the same exponents should be valid in the M-theory 
regime, where we consider the large- -/V limit with fixed 
9s ~ 9ym an d with the scaling \t — t'\ ~ N\t — r'\ of the 
light-cone time (in the string unit) , corresponding to the 
infinite momentum limit P + = N/g s — > oo. For possible 
interpretations of the predicted exponents from the 11 di- 
mensional viewpoint of the Matrix theory conjecture, we 
refer the reader to Ref. in addition to [9]. More de- 
tails of our analysis as well as results on other correlation 
functions shall be reported in a separate paper. 
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